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We develop a simple coarse-grained bead-spring polymer model exhibiting competing crystalliza- 
tion and glass transitions. For quench rates slower than the critical nucleation rate \T\ cr i t , systems 
exhibit a first-order crystallization transition below a critical temperature T = T crys t. Such systems 
form close-packed crystallites of FCC and/or HCP order, separated by domain walls, twin defects, 
and an amorphous interphase. The size of amorphous regions grows continuously as the quench rate 
\T\ increases, producing nearly amorphous structure for \T\ > |T| C rit- Our model exhibits many 
features observed in recent studies of crystallization of athermal polymer packings, but also critical 
differences arising from the softness of the pair interactions and the thermal nature of the phase 
transition. The model is considerably more computationally efficient than other recent crystalliz- 
able coarse-grained polymer models; while it sacrifices some features of real semicrystalline polymers 
(such as lamellar structure and chain disentanglement), we anticipate that it will serve as a useful 
model for studying generic features related to semicrystalline order in polymer solids. 

PACS numbers: 64.70.km,64.60.Cn,64.70.Dg,83.80.Ab 



I. INTRODUCTION 

A coarse-grained polymer model should include the 
minimal set of features necessary to capture the phys- 
ical phenomena of interest while remaining maximally 
computationally expedient. For example, the flexible 
Kremer-Grest (KG) bead-spring modelQ] is a minimal 
model including only chain connectivity, excluded vol- 
ume and van der Waals attractions. Despite this sim- 
plicity, it is able to capture the behavior of real polymers 
to an extraordinary degree, exhibiting features ranging 
from Rouse and entangled dynamics (i.e. reptation[2 ) 
in its molten state [U |3j to dynamical heterogeneity in 
its glass transition regime |4! to aging, rejuvenation, and 
strain hardening in its amorphous glassy state. pHT 

One limitation of the standard KG model is that it 
possesses an inherent length-scale competition; the equi- 
librium length £q of covalent bonds is significantly dif- 
ferent from the equilibrium separation tq for nonbonded 
monomers. This competition prevents formation of the 
semicrystalline order possessed by most real polymers. 
United atom models[Sl H] exhibit crvstallizationfTU] as 
well as glass formation. [llj and include the angular 
and dihedral interactions required to map to specific 
polymer chemistries, but are computationally expensive. 
In the opposite limit, the simplest models treat poly- 
mers as freely-jointed chains of tangent hard spheres 
with £q — tq and have recently illustrated the com- 
petition between athermal glass formation (jamming) 
and crystallization. [T^MTTj] The limitation of these latter, 
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highly-idealized models, of course, is that they are ather- 
mal, while fcgT is a critical parameter that profoundly 
affects polymer properties. It is desirable, therefore, to 
develop simple models which possess both the soft ex- 
cluded volume and van der Waals attractions necessary 
to capture thermal behavior (in particular, exhibiting a 
glass transition) and a local chain structure amenable to 
crystallization, e.g. £q = Vq. 

In this paper we develop and describe the basic proper- 
ties of such a polymer model. We will show that rapidly 
quenched systems remain largely amorphous down to 
T = while slowly quenched systems display a degree of 
crystalline order that increases with decreasing quench 
rate \T\. Consistent with results for athermal polymer 
packings[14-16 , our model forms close-packed crystal- 
lites of face centered cubic (FCC), hexagonal close packed 
(HCP), or mixed FCC/HCP order with varying degrees 
of stacking faults and five-fold-symmetric defects. While 
real semicrystalline polymers typically do not form close- 
packed crystals, we will show that our model captures 
generic features of polymer crystallization. 



II. MODEL AND METHODS 

Each polymer chain contains N — 50 coarse-grained 
beads, while each bead corresponds to 2-5 monomers. [1 
All beads have mass m and interact via the truncated and 
shifted Lennard- Jones potential Ui,j(r) = 4uo[(a/r) 12 — 
(a/r) 6 — (a/r c ) 12 + (a/r c ) 6 ], where a is monomer diam- 
eter, r c is the potential cutoff radius, and Ui,j(r) = 
for r > r c . The unit of time is r = \J ma? Juq and 
maps to time scales in the 10-100ps range; pQ we em- 
ploy a timestep St = r/300. We set r = £q = a 
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by choosing the pair interactions J/lj = 4uo[(<r/r) 12 — 
(cr/r) 6 - (cr/r c ) 12 + (cr/r c ) 6 ], with a = 2- 1 / 6 a and 
r c = 2 7 / 6 a, and using a stiff harmonic bond potential 
of form U bond (£) = (k b /2)(l - a) 2 , with k b = 600e. The 
energetic barrier to chain crossing is k b (\/2 — l) 2 ~ 100u , 
i.e. > 100/csT for the systems considered here. Systems 
consist of N c h = 500 chains and periodic boundary con- 
ditions are applied in all three directions. Intial melt 
states are generated with a monomer number density 
p = 1.0a -3 (packing fraction <fi = 7173/6). After thor- 
ough equilibration at fcsT = 1.2uq, systems are quenched 
to zero temperature at various rates \T\ while maintain- 
ing zero hydrostatic pressure using a Nose-Hoover baro- 
stat. Simulations are performed using LAMMPS.fTT] 
Throughout the rest of the paper, we will express tem- 
peratures in units of ksT/uo, quench rates in units of 
t _1 , distances in units of a and densities in units of a~ 3 . 

During the quenches we monitor several quantities in- 
cluding the potential energy per monomer U, pair cor- 
relation function g{r), packing fraction </>, and metrics 
of local structure including the Characteristic Crystallo- 
graphic Element (CCE) norm. 14- 16} [18] The later is a 
highly discriminating descriptor which quantifies the ori- 
cntational and radial similarity of a local environment to 
a given ordered structure in atomic and particulate sys- 
tems. The CCE norm is built around the defining set of 
crystallographic elements and the subset of distinct ele- 
ments of the corresponding point symmetry group that 
uniquely characterize the reference crystal structure. For 
example, the FCC crystal symmetry is mapped onto a set 
of four three- fold axes (roto- inversions of 27r/3), while 
the HCP is mapped onto a single six-fold symmetry axis 
(roto-inversion of 7r/3). A scan in the azimuthal and po- 
lar angles identifies the set of axes that minimize the CCE 
norm of a reference site (atom or particle) with respect to 
a crystal structure X. Details on the underlying mathe- 
matical formula and the algorithmic implementation can 
be found in Ref. [TH]- Once the CCE norm {ef) is calcu- 
lated for each site i a corresponding order parameter s x 
can be calculated which is practically equal to the frac- 
tion of sites with CCE norms below a pre-set threshold 
value (ef < e thrcs ). Results from the CCE-norm-based 
analysis with respect to FCC, HCP, and fivefold symme- 
tries are presented below. 

III. RESULTS 

Figure [l] illustrates the evolution of packing fraction 
4>{T) and potential energy U (T) at various quench rates. 
For slow quench rates, the data show clear signatures 
of a crystallization transition at T cryst ~ 0.56; <fi (U) ex- 
hibit upward (downward) jumps that indicate an increas- 
ingly first-order-like transition as \T\ decreases. As T ap- 
proaches zero, cj) approaches the maximal value for close 
packed crystals, <j) cp — n/y/18 = .7405. Indeed, at the 
slowest \T\, (f> exceeds <j) cp ; this is attributable to the soft- 
ness and long-range attractive tail of L/l j . In contrast to 



recent work on athermal systems [T2l 1 13] . our model does 
not "jam" at random close packing (<Prcp — .636: [15]) 
like atomic Lennard Jones systems [2D], it is an excellent 
crystal-former. For the fastest quench rate, 4> and U 
show no apparent crystallization transition, and a weak 
glass transition, as indicated by a smooth bend in U 
and (j), is observed at kgT = T g ~ 0.45ito. Thus, as in 
real semicrystalline polymers, according to the proposed 
model T cryst > T g . 
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FIG. 1: Volumetric and energetic measures of the crystalliza- 
tion transition. Solid blue, dotted green, and dashed red lines 
show data for \T\ = 1CT 6 , |T| = 1(T 5 , and \f\ = 10" 4 , respec- 
tively. Panel (a) illustrates the packing fraction <j> and panel 
(b) illustrates the potential energy per monomer U. For the 
slowest quench rates, both data sets indicate T crys t — 0.56. 

Figure [2] shows the evolution of the pair correlation 
function g(r) with T at the slowest and fastest \T\. Re- 
sults are shown for temperatures well above the melting 
point, slightly below T cryst , and zero. Above the melting 
point, systems have amorphous (melt-like) structure as 
expected. For slow quenches, just below T cryst , peaks in 
the correlation function form corresponding to the ap- 
pearance and growth of close packed order. At zero tem- 
perature, clear peaks at the characteristic second and 
third nearest neighbor distances for close-packed crystals, 
T2n = \/2 and r^ n — V% have developed|21j: the system 
also retains some amorphous character as indicated by 
the large width of these peaks. In sharp contrast, for 
fast quenches |T| = 10 -4 , systems at the same tempera- 
tures remain predominantly amorphous; g(r) maintains 
liquid-like structure down to T = 0. 

We now turn to a detailed examination of local envi- 
ronment around each monomer and to the identification 
of the crystalline structure (or the lack thereof) at vari- 
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noncrystalline portion of samples consists of stack- faulted 
domain walls and an amorphous "interphase" analogous 
to that found in real semicrystallinc polymers. 



FIG. 2: Pair correlation functions at characteristic tempera- 
tures and different quench rates. Red, purple and blue lines 
respectively show data for knT/uo = 1.0, 0.5 and 0. Panel 
(a) shows data for \T\ — 10 -4 while panel (b) shows data for 
If I = 10" 6 



ous T, mainly in the regime of T cryst and at T = 0. As 
described in the methods section, we have implemented 
the CCE norm and compared against the HCP, FCC, and 
fivefold structures. The highly discriminating nature of 
the CCE norm is demonstrated in Figure [3] where parity 
plots|18) for the HCP and FCC CCE norms are shown for 
|T| = 10~ 6 at T — 1.0, where the system is amorphous, 
and at T = where it becomes predominantly ordered. 
Sites with HCP- (or FCC-) CCE norms with values lower 
than e thros = 0.20 are characterized as HCP-like (or FCC- 
like). By construction, a monomer with high HCP simi- 
larity (i.e. low HCP-CCE norm) possesses low FCC sim- 
ilarity (high value of FCC-CCE norm) and vice versa. 
Thus, for any system configuration we can reliably iden- 
tify the local environment around each monomer with 
respect to HCP, FCC and fivefold symmetries. The fig- 
ure illustrates that in the liquid state at T = 1.0, es- 
sentially no monomers have local FCC or HCP order, 
while at zero temperature, a very high fraction of sites 
(~ 65%) possess either FCC or HCP order, with compa- 
rable probability. The remaining 35% of sites have either 
five-fold or other local structure, as indicated by the re- 
gion (e FCC , e HCP ) > (0.2, 0.2). In parallel, the vacancy of 
the region defined by (e FCC ,e HCP ) < (0.2,0.2), as seen 
in both panels, highlights the discriminating character of 
the CCE-norm descriptor. As we will show below, this 
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FIG. 3: Parity plot of the e FCC versus e HCP CCE-based norms 
over all monomers for \T\ — 10 -6 at T — and (inset) 
T = 1.0. Colors indicate the (log-scale) site ordering probabil- 
ity density P(e HCP , e FCC ). Horizontal and vertical gray lines 
indicate the threshold value of the CCE norm (e thrcs = 0.20). 

Next we present results, based on the CCE analysis, 
on the evolution of local ordering during cooling. Fig- 
ure [4] shows the fraction of sites with: (a) close-packed 
(FCC or HCP) order, (b) fivefold similarity, and (c) nei- 
ther fivefold local symmetry nor close-packed order, as a 
function of T for the various quench rates. In all cases, 
close-packed ordering grows continuously as T decreases, 
with the transition at T = T cryst becoming increas- 
ingly first-order-like with decreasing |T|. As T continues 
to decrease, the fraction of close-packed sites continues 
to increase, indicating an effective "annealing" process 
wherein defects are removed. Throughout this process, 
the fraction of sites with FCC order is comparable to but 
exceeds the fraction of sites with HCP order, especially 
at the slowest \T\. This is expected, since while the free 
energy difference between FCC and HCP phases is very 
small[22], crystal-growth kinetics favor FCC crystallite 
formation. 

The fraction of sites with close-packed order, f cp , in- 
creases sharply with decreasing quench rate for T slightly 
below T cryst , and continues to increase as T decreases 
to zero. For example, at T = 0, f cp is only 10% for 
|f I = 10~ 4 , but 58% for \f\ = 10~ 5 and 65% for 
|T| = 10~ 6 . At the intermediate quench rate, the jumps 
in U, (j), and f cp all exhibit a "delay" to T ~ .52, indi- 
cating that the critical nucleation rate |T| cr j4 ~ 10~ 5 . 
Quench-rate-dependent differences in crystal structure 
for T — will be examined in more detail below. 
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FIG. 4: Measures of crystalline order vs. T at various quench 
rates. Solid blue, dotted green, and dashed red lines show 
data for |T| = 10~ 6 , |T| = 10~ 5 , and \f\ = 1CT 4 , respectively. 
Panel (a) : fraction of sites f cp with close-packed order, Panel 

(b) : fraction of sites f 5 f with fivefold local symmetry, Panel 

(c) : fraction of sites foth with other local structure. 



Fivefold local symmetry is well-known to inhibit crys- 
tallization and promote amorphous structure. [24 29 For 
all quench rates, it is clear that as density (<j>) increases 
for T > T crystl so does the population of sites with 
fivefold symmetry. This trend is in perfect agreement 
with past findings from simulations on monomeric hard 
spheres of uniform size|29j where fivefold probability in 
amorphous packings increases as the system becomes 
denser. The physical trend changes drastically as tem- 
perature reaches and drops below T cryst . For the fastest 
quench rate, fivefold sites continue to grow linearly and 
the system remains amorphous with only a small frac- 
tion of ordered sites. In contrast, the population of five- 
fold sites drops significantly for \T\ = 1CP 5 , and remains 
nearly constant for \T\ = 10~ 6 . This finding clearly 
points towards a structural competition between close- 



packed ordering and fivefold symmetry (e.g. twin de- 
fects), a physical trend that has also been observed in 
a wide range of athermal hard-sphere packings. [Z5M29"] ■ 

Many sites lack either close-packed order or five-fold 
similarity; panel (c) shows the fraction of such sites, 
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foth 

indicating that for the CCE-norm structure-identification 
procedure described above, f th is a good discriminant of 
liquid-like order. For |T| > |T| cr i t , systems retain amor- 
phous structure down to T — 0, consistent with the g(r) 
data shown in Fig. [ij For \T\ < \T\ cr it, foth shows a 
first-order-like (downward) jump at T = T cryst , and con- 
tinues to decrease with decreasing T during "annealing" , 
but remains significant down to T = 0. Our model is 
therefore well-suited to producing the semicrystallinity 
observed in real polymers. 




FIG. 5: Snapshot of system quenched at \T\ = 10 -6 , at 
T — 0. (Top panel) HCP-ordered sites are shown in blue, 
FCC-ordered sites in red, and fivefold sites in green. (Bottom 
panel) "Other" sites are shown in green. Image created with 
the VMD software. [30] 



Visualization of systems prepared at various \T\ pro- 
vides considerable insight into the semicrystalline mor- 
phologies formed by our model. Figure [5] shows the end 
state (T = 0) of the slowest quench. Grain-like HCP 
and FCC domains are clearly visible, and the fivefold- 
symmetric sites often correspond to twin defects - a struc- 
ture similar to that found in polycrystalline metallic or 
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colloidal systems (see e.g Ref. [3T]) and model packings of 
monomeric hard spheres. [28] [29] The ordered structures, 
as established here, show reduced tendency to layer for- 
mation (randomly stacked hexagonal close packing) com- 
pared to that found for hard-sphere chains in Refs. fPfl - 
fT6] , presumably because the larger system sizes employed 
here reduce the influence of the periodic boundaries or 
because in the athermal systems a strict tangency condi- 
tion is applied with respect to bond lengths. Gaps cor- 
responding to the "other" sites are clearly visible in the 
snapshot. In the bottom panel, we illustrate these sites 
for the same system. Careful visual inspection shows 
that these regions possess nearly close-packed structure 
and correspond to stack-faulted domain walls. Since the 
crystallite domain size is considerably smaller than our 
simulation cells, these domain walls form a percolating 
structure. 




FIG. 6: Snapshot of system quenched at \ f\ = 1CT 4 , at T = 0. 
HCP-ordered sites are shown in blue, FCC-ordered sites in 
red, and fivefold sites in green. 

Faster quench rates produce reduced crystalline and 
greater amorphous order. Figure [6] shows the (T = 0) 
end state of a system quenched at \T\ = 10~ 4 . HCP 
and FCC crystallites are present, but are far smaller 
and fewer. Furthermore, the crystallite/grain-boundary 
structure produced for |T| < |T| cr it is absent for fast 
quench rates. In parallel, the number of fivefold sites 
is much greater than for \T\ < \T\ crit , and these sites, 
rather than corresponding to twin defects, are apparently 
arranged randomly. Visual inspection of the "other" 
sites for this quench rate shows that they are much less 
ordered than those for the lowest quench rate, in ef- 
fect corresponding to an amorphous interphase like that 
found in real semicrystalline polymers |32j . Thus our 
model is able to produce large crystallites with domain 
walls for |T| < |T| cr i t and a predominantly amorphous 
structure with small crystallites for \T\ > \T\ cr i t . For 
\T\ = 10~ 5 ~ \T\ crit , results are intermediate between 
these two limiting cases, with a tendency towards the 
crystalline state. 

Finally, significant information can be gained on crys- 
tal nucleation and growth can be obtained by visual ex- 
amination of systems at T ~ T cryst . Figure [7] shows 




FIG. 7: Crystal nucleation at |T| = 1(T6. (Top) Formation 
of a preliminary nucleus at T = .562. (Bottom) growth of nu- 
cleus: T = .561. The systems have f cp = 0.019 and 0.14, re- 
spectively. HCP-ordered sites are shown in blue, FCC-ordered 
sites in red, and fivefold sites in green. Image created with 
the VMD software [3D]. 

a series of snapshots from the \T\ = quench. At 

T = 0.562 the first trace of crystal aggregates can be 
seen in the form of a "baby" nucleus which consists of 
similar amounts of HCP and FCC sites. As the system is 
still amorphous, the number of sites with fivefold symme- 
try is comparable to the fraction of sites with either HCP 
or FCC similarity. However, by T = 0.561 the number of 
ordered sites present in the system has greatly increased 
and the first large crystal seed (critical nucleus) is extant 
consisting again of roughly equal amounts of HCP- and 
FCC-like sites. As T continues to drop, this nucleus con- 
tinues to grow and expand until it fills most of the system 
as illustrated in Fig. [5] 

IV. DISCUSSION AND CONCLUSIONS 

We have developed a simple, computationally effi- 
cient bead-spring model exhibiting competing crystal- 
lization and glass transitions. Crystallization is pro- 
moted by removing the length scale competition present 
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in the Kremer-Grest bead-spring model. pQ At quench 
rates faster than the characteristic nucleation rate |r| cr jt, 
systems remain predominantly amorphous with a large 
number of fivefold-symmetric sites, typifying quenched- 
disorder vitrification. At quench rates slower than |T| cr .j t , 
the system exhibits a first-order-like phase transition 
wherein crystal nuclei form in coexistence with an amor- 
phous phase. Since our model employs flexible chains and 
possesses a covalent bond length that is equal to the bead 
diameter, the crystals formed by for slow quench rates 
exhibit mixed HCP/FCC close-packed order. While an- 
gular interactions and competing length scales produce 
different crystalline structures in united-atom-model sim- 
ulations of chemicaly specific polymers, the close-packed 
order formed by our systems may be viewed as a conse- 
quence of the level of coarse-graining; as in the Kremer- 
Grest model, one bead corresponds to several monomers. 
We expect that this model should be useful in studies 
of how semicrystalline order affects phenomena such as 
aging, dynamical heterogeneity, and the nonlinear me- 
chanics of solid polymers. 

Our model is comparable to but simpler and compu- 
tationally "cheaper" than a recent bead-spring model of 
polyvinyl alcohol (PVA) [33J IS] which has also been used 
to study generic features of polymer crystallization. The 
model of Refs. [331 E3] generates more realistic features 



of crystallization such as formation of lamellae and chain 
disentanglement during the crystallization process, but 
the lamellar equilibrium crystal structure generated by 
the stiff angular potential employed in the model makes 
its critical quench rate for crystal formation about 1.5 
orders of magnitude (i.e. ~ 10~ 65 ) lower than that re- 
ported here. Here we have sacrificed some chemical re- 
alism to achieve a minimal model of thermal polymer 
crystallization. Future work will develop greater chem- 
ical realism and the ability to form lamellar structures, 
without sacrificing computational efficiency, by integrat- 
ing generic angular potentials. 
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